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ABSTRACT 


t 

A  numerical  procedure  for  the  solution  of  the  nonlinear  equations  govern¬ 
ing  the  large  axisymmetric  deflections  of  thin  shells  of  revolution  is  pre¬ 
sented  and  applied  both  to  the  complete  equations  due  to  Reissner  and  to  the 
simpler  equations  to  vhich  these  reduce  for  small-finite  angle  changes. 

Global  solutions  extending  into  the  postbuckled  range  are  shown  to  be  consid¬ 
erably  more  complicated  than  expected.  The  character  of  the  global  solution 
is  also  shown  to  be  quite  sensitive  to  boundary  conditions  imposed.  A  compar¬ 
ison  of  the  results  obtained  from  the  complete  equations  and  the  small-finite 
deflection  equations  reveals  a  very  close  agreement  through  the  entire  load- 
deflection  history. 
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NOMENCLATURE 

slope  of  meridian  on  deformed,  undeformed  middle  surface 
shell  thickness 

horizontal  and  vertical  coordinates  of  undeformed  middle  surface 

(r;)2  ♦  (z^)2,  where  primes  refer  to  differentiation  with 
respect  to  the  independent  variable,  £. 

Poisson's  ratio.  Young's  modulus 
l/(Eh) 

Eh3/12(l-v2)) 

vertical,  horizontal  stress  resultants 

r  H  *  stress  function 
o 

Uo  -  4)  ■  rotation  of  meridional  tangent 

uniform  pressure,  psi 
nondimensionalized  pressure 

horizontal,  vertical  components  of  uniform  pressure 

concentrated  load  at  apex  of  spherical  shell 
aP/(2»Eh3)  *  nondimensional  load  parameter 
meridional,  circumferential  stress  resultants 
meridional,  circumferential  stress  couples 
horizontal,  vertical  deflections 
vertical  deflection  at  apex  of  sphere 
mesh  spacing 

radius  of  spherical  shell 
half-angle  opening  of  spherical  shell 

12(l-v2) 

m2C2  a/h  «  nondimensional  geometric  parameter  for  shallow 
spherical  shell. 


iii 


INTRODUCTION 


The  objective  of  this  paper  is  the  description  and  application  of  a 
numerical  procedure  for  the  solution  of  the  nonlinear  differential  equations 
governing  the  finite  axisymmetric  deformation  of  thin  shells  of  revolution  as 
presented  by  Reissner.1  The  procedure  is  applied  both  to  the  complete  equa¬ 
tions  (I  and  II  of  Reference  l),  valid  for  arbitrarily  large  deflections  con¬ 
sistent  vith  small  strains,  as  veil  as  to  the  simpler  set  of  equations  (III 
and  IV  of  Reference  l)  valid  for  small-finite  angle  changes.  The  latter  set 
of  equations  (or  their  equivalent)  has  formed  the  basis  for  many  previous 
investigations  in  this  area.  Thurston7  has  studied  the  complete  equations 
but  does  not  report  on  a  comparison  of  the  results  of  the  two  sets  of  equa¬ 
tions.  The  present  study  consists  essentially  of  an  application  of  Newton' s 
method8  to  obtain  a  system  of  linear  correctional  equations  for  an  initial 
approximate  solution  and  the  employment  of  a  Gaussian  elimination  procedure9 
for  the  solution  of  the  finite  difference  equivalent  of  this  linear  system. 
Archer,4  Wilson  and  Spier,10  and  others  have  previously  employed  a  Gaussian 
elimination  procedure  in  this  problem  area,  while  Thurston5  has  utilized 
Newton's  method  for  the  case  of  uniform  pressure  on  a  clamped  spherical  cap. 
The  specific  advantage  provided  by  the  procedure  discussed  in  this  paper  is 
that  it  more  completely  characterizes  the  solutions.  In  particular,  it  per¬ 
mits  development  of  global  solutions  which  are  shown  to  be  continuous  from 
the  small  deflections  encountered  in  the  prebuckled  state  through  the  large 
deflections  of  the  postbuckled  state.  The  continuous  character  of  these 
solutions  is  not  only  interesting  from  a  theoretical  point  of  view  but  is 
also  useful  in  permitting  one  to  develop  a  complete  load-deflection  history 
for  a  new  geometry  and  loading  condition  without  having  a  priori  estimates  of 
the  solution  in  either  the  prebuckled  or  the  postbuckled  state. 

To  illustrate  this  latter  point,  we  consider  the  problems  of  a  clamped 
spherical  shell  under  uniform  pressure  and  that  of  a  spherical  shell  under  a 
concentrated  load  at  the  apex,  with  a  clamped  edge  and  vith  an  unrestrained 
edge.  New  results  are  presented  for  these  problems  in  the  so-called  post¬ 
buckled  range. 
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BASIC  EQUATIONS 


Reissner1  has  given  the  equations  governing  the  finite  axi  symmetric 
deflections  of  thin  shells  of  revolution  in  the  form 

♦  ^o/F0)U-^0)'  -  (oQ/ro)2(cos  *)(sin  +  -  sin  #o) 

♦  v  (o Jr  )[(cos  +  -  cos  ♦  )♦'  ♦  (D'/DHsin  +  -  sin  ♦  )] 

*  a*/(r  D)  sin  +  -  r  V  cos  ♦]  (1) 

♦  (Gq/Go)i>'  -  [(a0  cos  +/rQ)2  -  v(ao/ro)(*'  sin  ♦  ♦  C'  cos 

♦  (a*C/r  )(cos  +  -  cos  ♦  )  ♦  vo  (sin  $)(r  V)'/r 

O  O  0  0  o  o 

♦  [(oQ/ro)2  cos  ♦  sin  ♦  ♦  v(ao/ro)U'  cos  ♦  -  c'  sin  ♦/C)]  rQV 

-  (oo/ro)(ro  V  -  Cv(«o/ro>2  co*  ♦  -  <«oC')/(roC)]  ro  pH  (2) 

wher*  F0  -  (r0D/o0).  O0  -  r^C), 

where  primes  denote  differentiation  with  respect  to  the  independent  variable 
£•  and  where  subscript  zero  refers  to  the  value  of  the  subscripted  quantity- 
before  deformation  (see  Figure  1). 

The  appropriate  stress  and  moment  resultants  and  displacements  are 
expressed  in  terms  of  ♦  and  i|*  by  the  relations 

roV  *  -  f rQPy  aod* 

r  K.  «  ji  cos  ♦  +  r  V  sin  a 

O  t,  o 

r  Q  *  -  ip  sin  4  +  r  V  cos  $ 
o  o 

“oNe  ■  *'  +  ro°opH  t3) 


2 


s  "d  [U-$0)'/ao  +  v(sin  <f>  -  sin  $Q)/rJ 

Mq  s  "D  [(sin  ♦  -  sin  ♦qJ/Fq  ♦  v(<J>-4>o)  ’/aj 
u  *  r  (N  -  vN  )/C 

o  O  ^ 

w  *  /[(sin  4  -  sin  4  )  ♦  sin  4  (II  -  vN  )/cl  a  d£. 

For  small- finite  deflections  (B2  =  UQ-$)2  «1)»  Equations  1  and  2 
simplify  to 


Fo8"  +  Fi6'  *  F2B  +  V  *  ri 


0  (j/11  ♦  G  i|> 1  +  C  ♦  G  B  *  r 
o  x  2  3  2 


(M 

(5) 


where 


F  m  F  1 
1  o 


*  Go’ 


F  «  -  (r')2  D/(r  a  )  ♦  v(r  D/a  )' 

2  O  O  O  O  O 


(r’)2/(r  a  C)  -  v(r'/o  C) * 

O  O  O  O  O 


F,  *  a  sin  $ 

3  o  o 


°3  *-F3 


r,  *  [(3r'i'D)/(?r  a  )  -  v(i'D/2o  )']S2 

1  u  O  O  oo  o  0J 

+  a  r  V  cos  $  ♦  a  sHp  cos  $  +  r  V  sin  $  ] 

OO  OO1-  OO  O'1 

r,  *  B2Vo)/(roaoC)  *  v(z^/aoC)']Blp  ♦  vz^B'<j»/(aoC) 


o  o 


(6) 


-  a  82cos  i  /2  +  [rV/(r  a  C)  ♦  v(z'/a  C) ']  (r  V) 
o  o  oooo  oo  o 

♦  v(z’/a  C)(r  V)'  ♦  [(z'2  -  r’2)/(o  r  C)  -  v(r '/a  C) ']  (r  V)B 

Q  Q  n  »  ft  n  O  n  A  n  —  A 


o  o 


o  o 


-  vr'(Br  v)7(a  c)  -  (r  2pH)’/C  -  [v(r’  +  Bz')  -  rc'/c](rpJC) 

OO  O  On  wO  O  O  OH 


The  corresponding  simplification  of  Equations  3  will  not  be  reproduced  here 
in  the  interest  of  brevity. 
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METHOO  OF  SOLUTION 


As  mentioned,  the  numerical  technique  employed  in  this  report  combines  the 
standard  Gaussian  elimination  procedure  with  the  classical  Newton  iteration 
method.  The  former  substantially  reduces  the  demands  made  upon  a  digital 
computer  with  respect  to  both  memory  and  speed.  The  latter  improves  conver¬ 
gence  significantly  an  I, when  coupled  with  a  fairly  direct  method  for  obtain¬ 
ing  starting  solutions ,  permits  the  development  of  solutions  in  the  so-called 
postbuckled  range. 

Details  of  the  procedure  are  perhaps  most  compactly  presented  in  terms 
of  the  small- finite  deflection  equations.  Extension  to  the  more  general 
equations  follows  easily.  The  essence  of  the  Newtonian  iterative  scheme  as 
outlined  in  Reference  8  consists  in  replacing  the  nonlinear  differential  equa¬ 
tions  by  a  sequence  of  linear  differential  equations.  Specifically,  a  cor¬ 
rection  ( 66j ,  fiiJ'j)  to  the  approximate  solution  (Bj,  i|/j)  is  sought  according 
to  the  prescription 


6  ■  6j  +  46j 

(7) 

*  *  *j  * 


where  (B,i|»)  is  the  actual  solution.  Inserting  (7)  into  (1*)  and  (5)  and 
omitting  terms  in  T  ,  T2  which  are  negligible,  the  result  may  be  written 

F  6B"  ♦  F.6B!  ♦  (F  -  a  i|>  cos  ♦  )  6B,  ♦  (F  -  a  B.  cos  $  )  6i|>, 
°J‘J  20J  oj  3  0,1  o.l 

.  r*  •  r,  -  {FoBj  ♦  FlB;  *  F2Bj  ♦  f3*j} 

G  ♦  G, 6^1  +  0  .  ♦  (G  ♦  a  B,  cos  4  )  56. 

°  J  1  J  2  J  3oJ  oj 

1  r2*  *  T2  '  <°o*J  +  Vj  +  Vj  +  G3  6jK 


(8) 


In  Equations  8  we  have  omitted  nonlinear  terms  in  6 Bj  and  ,  since  for  am 

assumed  solution  (B.,  4».)  reasonably  close  to  the  correct  solution  these  te: 

J  J 

are  small  compared  to  the  linear  terms.  An  iterative  procedure  is  now 


U 


adopted  in  which  Equations  8  are  solved  for  ( *  6^)  and  a  new  approximate 
solution  'J'j+j)  *8  o^a^ned  from 

V.  *  8j  +  <bj 

(9) 

Vi  ’  *J +  sv 

If  the  correction  terms  (68.,  6i|/  )  approach  zero  as  the  number  of  iterations 

J  J 

increase,  then  the  approximate  solution  obtained  by  this  iterative  process 
approaches  an  exact  solution  of  the  original  equations.  For  the  solution  of 
Equations  8,  we  replace  derivatives  by  simple  central  differences  over  a  mesh 
of  n  points  with  spacing  k  and  obtain 


Ai<Ti-.  ♦  Bi  STi ♦  ci  <Ti+1  ■  Di 


(10) 


Jl  *  2|3|  •  •  •  t  n-1 


where,  in  the  compact  notation  employed  by  Archer,4 


(11) 


F  -  kF,/2  0 

o  1 


oo  -  »,/*/ 


•  c. 


F  ♦  kF,/2  0 

o  I 


°0  +  k0l/2/l 


-2Fo  +  k2(F2  -  «o(eos  *0)*)  k2(Fs  -  o0(co»*0)B) 


k2(G  ♦  a  (cos  ♦  )8) 
30  o 


-2G  +  k2G_ 

O  2 


ru 


(12) 


If  the  original  boundary  conditions  are  formulated  as 


B1T1  *  C1T2  *  Dl 


(13) 


AT  .  ♦  B  T  ■  D 
n  n-1  n  n  n 
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then  the  boundary  conditions  for  the  modified  problem  become 

B|«l  +  C1«T2  *  DJ  -  Vl  ’  C!T2  *  D.#  *  0 


A  «T  ,  +  B  6T  ■  D  -  AT  -  B  T  ■  D  *0. 
n  n-1  n  n  n  n  n-i  n  n  n 


(lk> 


The  process  of  Gaussian  elimination  very  efficiently  inverts  this  system 
of  equations.  This  process  may  be  compactly  summarized  by  the  relations4 


M.  ■  Bi 


S,  •  WJ'd, 


R,  -  W71C, 


Hi  *  Bi  -  AiRi-l 


Ri  *  Milci 


W  *  B  -  A  R  , 
n  n  n  n-1 


si*wI1(Di  -Vi.,) 


i  *  2 13| • • • •  n-1 


s  *  W_l(D  -AS  ,  ) 
n  n  n  n  n-1 


(15) 


T  *  S  T,  *  S.  —  R.T  .  i  *  n-1,. . . .  1. 

n  n  i  i  i  i-1 


A  detailed  description  of  the  selection  of  starting  values  and  of  the  pro¬ 
cedure  for  obtaining  solutions  in  the  postbuckled  regions  is  best  given  with 
reference  to  a  specific  example  and  is  postponed  until  the  next  section. 

Reissner's  complete  equations  may  be  treated  in  essentially  the  same 
fashion,  vith  terms  such  as  cos  d  replaced  by 

cos  ♦  =  cos  (d  -  6)  =  cos  ($  -6.-66.) 

o  oil 

«  cos  (♦  -  6.)  +  66,  sin  U  -  8.)  (l6) 

O  1  1  o  1 

and 

sin  ♦  =  sin  (d  -  6.)  -  66.  cos  (d  -  6.). 

O  1  1  O  i 


The  finite  difference  simulation  of  the  complete  equations  then  has  the  same 
form  as  that  for  the  small- finite  deflection  equations.  The  modifications  of 
the  definitions  for  FK»  are  not  presented  here  in  the  interest  of  brevity, 
and,  also,  in  view  of  the  outcome  of  the  comparison  to  be  made  between  the 
results  of  the  two  sets  of  equations. 
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NUMERICAL  RESULTS 


a.  Concentrated  Load  on  Spherics'  Shell 

Consider  first  the  case  of  a  spherical  segment  vith  a  concentrated  load  P 
at  the  apex  and  an  unrestrained  edge.  For  this  problem  the  equations  of  the 
middle  surface  of  the  shell  are  taken  as 


while 


r  *  a  sin  z  *  -a  cos  £ 

o  ’  o 


PH  *  Pv  *  °»  (rQV)  *  P/2*. 


(IT) 

(18) 


Convenient  nondimensional  geometric  and  load  parameters  are  defined  as 

X2*  m2£2  a/h  and  P  *  aP/(2irEh3). 
o 


(19) 


We  observe  that  although  it  is  not  necessary  to  confine  attention  to  shallow 
shells,  nonetheless,  due  to  the  nature  of  the  load,  this  is  the  most  interest¬ 
ing  area  of  application.  If  deep  shells  are  being  considered,  the  definition 
of  A2  involves  sin2(£Q)  rather  than  £2  .  Appropriate  boundary  conditions  are 


8  *  <>  *  0  at  £  *  0 

*  0  at  £  * 


(20) 


With  this  information  one  may  proceed  as  follows:  For  a  given  value  of 

« 

X  and  P  the  linear  solution  is  first  obtained  by  suppressing  nonlinear  terms 

.  * 

in  Tj  and  r2.  (In  this  connection,  P  is  initially  chosen  small  enough  that 

the  linear  solution  is  reasonably  appropriate).  The  linear  solution  is  used 

,  * 

as  an  initial  estimate  of  the  nonlinear  solution  (8  ,  )  for  the  same  P  , and 

oo 

Equations  8  are  solved  for  the  corrections  (6B0,  6^Q)  to  this  initial  esti¬ 
mate.  The  corrected  set  of  solution  values  (Bj,  i|>  )  is  used  to  obtain  a  new 
set  of  corrections  (68j,  6^),  and  the  process  is  repeated  until  the  solution 
converges  to  a  specified  degree  of  accuracy.  In  this  connection,  it  should 
be  noted  that  as  a  convergence  criterion  we  employed  the  requirement  that  for 
the  ftph.  iteration  and  for  each  point  i  on  the  mesh 
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* 

When  satisfactory  convergence  has  been  achieved  for  a  given  value  of  P  , 

* 

stresses  and  displacements  along  a  shell  meridian  may  be  calculated.  Then  P 

is  incremented  and  the  previously  converged  solution  for  the  lower  value  of 

#  * 

P  is  used  to  start  the  iteration  for  the  higher  value  of  P  .  In  this  manner 

curves  of  load  versus  deflection  are  obtained  for  given  values  of  X.  It  is  a 
relatively  direct  matter  to  move  along  a  branch  of  the  load-deflection  curve 
until  a  local  maximum  (or  minimum)  is  reached.  At  such  a  point,  the  following 
simple  procedure  was  found  to  be  adequate  for  obtaining  starting  solutions  on 
the  neighboring  (continuous)  branch.  For  a  value  of  P‘  slightly  below  the 
maximum  (or  above  the  minimum),  take  as  a  starting  set  of  B  and  ii  a  multiple 

/  X  * 

(9j3,  02<J/)  of  the  solution  found  for  P  on  the  preceding  branch.  A  very 

small  amount  of  numerical  experimentation  is  sufficient  to  find  constant  values 
of  0  and  0  which  produce  a  convergent  solution  on  the  new  branch.  These 

*  n 

values  of  0j  and  02  depend  on  A2  and  P  ,  but  generally,  1.0  <  (0j,  ©2)  <  3.0 
for  branches  moving  to  the  right  (larger  deflections),  while  0.50  <  (0j,  @2)  < 
1.0  for  branches  moving  to  the  left  (smaller  deflections).  It  is  now  a  simple 
matter  to  move  along  the  new  branch  until  another  local  minimum  (or  maximum) 
is  found. 

For  the  specific  problem  under  discussion,  the  load  deflection  curves 
for  small  A  (say  A2  5  lU)  are  found  to  be  monotonically  increasing.  As  A 
increases,  local  maxima  and  minima  emerge  but  the  curves  still  retain  the 
generic  shape  frequently  associated  with  (or  postulated  for)  the  mechanism  of 
buckling  (sec  Figure  2).  Our  objective  in  this  report  is  not  so  much  the 
prediction  of  buckling  loads  for  this  problem  (this  was  done  in  Reference  11, 
where  excellent  agreement  with  the  experimental  results  of  Evan-Iwanowski ,  et 
al,  was  demonstrated) , but  is  concerned  rather  with  the  observation  that  for 
higher  A  values  the  character  of  the  load-deflection  curve  changes  consider¬ 
ably.  For  example,  in  Figures  3a,  b,  and  c,  we  present  numerical  results  for 
A?  =  8l,  and  lUU. 


To  emphasize  the  dependence  of  the  solution  upon  boundary  conditions  we 
present  in  Figure  k  the  load-deflection  curve  obtained  for  the  same  problem, 
with  X2  *  lUU  but  with  a  clamped  edge  and,  therefore,  the  boundary  conditions 


8  « 

(j,  . 

0 

at 

ii 

o 

8  * 

u  * 

0 

at 

It 

o 

The  difference  in  behavior  of  the  solution  produced  by  this  change  in  the 
boundary  condition  is  rather  surprising  in  view  of  the  loading.  We  observe 
that  there  is  also  experimental  evidence  of  a  marked  difference  in  the 
behavior  of  clamped  versus  unrestrained  spherical  caps  under  concentrated 
loads.  Evan-Iwanowski, et  al,12  concluded  on  the  basis  of  detailed  experimental 
studies  that  such  a  clamped  shell  exhibits  the  load  deformation  pattern  shown 
in  Figure  U  and,  therefore,  does  not  buckle. 

Finally,  we  observe  that  the  numerical  results  obtained  by  using  the 
complete  equations  differed  only  very  slightly  from  those  obtained  from  the 
small- finite  equations.  The  overall  character  of  the  two  sets  of  results  is 
the  same,  i.e.,  the  load-deflection  curves  have  the  same  continuous  character 
and  the  same  number  of  branches.  Numerically,  the  typical  deviation  along  a 
branch  of  the  curves  is  of  the  order  of  one  percent,  even  for  the  final  post- 
buckled  branch  where  maximum  discrepancy  is  to  be  expected.  The  greatest 
deviations  occur  near  the  maxima  or  minima  of  a  branch,  but  even  there  the 
discrepancy  is  so  slight  that  the  difference  between  the  load-deflection 
curves  is  barely  discernible  on  a  graph  of  reasonable  scale. 

b.  Uniform  Pressure  on  Spherical  Cup 

Turning  next  to  the  case  of  uniform  pressure  p  on  a  clamped  spherical 
shell,  we  have  for  loading  conditions 

2 

PH  *  p  sin  £•  py  *  -p  cos  K  rQV  =  rQ  p/2  (23) 

where 

rQ  =  a  sin  £  zq  =  -a  cos  £.  (2*0 

Nondimensional  geometric  and  load  parameters  are 
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p  *  -m2a2p/(^Eh2) . 


(25) 


A2  *  m2S2  a/h  and 
o 

We  shall  also  make  use  of  the  "average  deflection"  parameter  p,  defined  as 

p  *  m2w/h 


where  R 

w  ■  (2/R2)  /  °  r  w  dr 
o 

o 


and 

R  * 

o 

Boundary  conditions  for  the  clamped 

8  *  *  0 

6  *  u  =  0 


a  sin  €  . 
o 

spherical  shell  are 
at  C  a  0 

at  £  *  CQ. 


(26) 


(27) 


Typical  load-deflection  curves  obtained  for  this  problem  are  illustrated  in 
Figures  5  through  8.  These  curves  bring  out  a  gradual  transition  from  the 
simpler  behavior  for  A  *  5  through  a  more  elaborate  behavior  for  A  *  8,  a 
return  to  a  relatively  simple  pattern  for  A  *  12,  and  finally,  the  emergence 
of  the  more  intricate  behavior  again  at  A  *  20. 

As  in  the  case  of  the  concentrated  load,  use  of  the  complete  equations 
produced  no  essential  change  in  the  qualitative  behavior  of  the  solution, 
viz.,  the  load-deflection  curves  are  continuous  and  have  the  same  number  of 
branches  as  do  the  solutions  of  the  small- finite  equations.  Quantitatively, 
the  agreement  is  again  very  good,  with  the  greatest  discrepancies  being  of  the 
order  of  a  few  percent  and  occurring  at  the  maxima  or  minima  of  the  curves. 

The  agreement  persists  even  into  the  final  branch  of  the  curves,  which  cor¬ 
responds  to  a  nearly  inverted  shallow  shell. 


DISCUSSION  OF  RESULTS 


The  present  numerical  results  demonstrate  the  existence  of  a  large  num¬ 
ber  of  distinct  equilibrium  positions  for  a  given  load  and  for  certain  ranges 
of  the  geometric  parameter.  The  existence  of  these  multiple  solutions  is  also 
shown  to  be  dependent  to  a  great  extent  upon  the  boundary  conditions  but  to 
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be  independent  of  the  employment  of  the  complete  equations  rather  than  the 
small- finite  deflection  equations.  Recently  (during  the  preparation  of  this 
paper,  in  fact),  it  came  to  the  writer's  attention  that  Anselone,  Bueckner, 
Johnson,  and  Moore,13'14  have  studied  the  large  deflections  of  a  clamped 
shallow  spherical  cap  under  uniform  pressure  and,  employing  a  completely  dif¬ 
ferent  technique  from  the  one  described  here,  have  obtained  results  which 
exhibit  the  same  continuous  character  and  which  agree  quite  well  numerically 
with  those  presented  here.  We  note  also  that  Thurstons  and  Keller  and  Reiss6 
have  pointed  out  the  possibility  of  multiple  equilibrium  configurations  for 
pressure  loading  on  a-  clamped  spherical  cap,  although  their  results  do  not 
bring  out  the  continuity  of  the  solutions. 

The  increasing  complexity  in  the  structure  of  the  solutions  as  X  increases 
brings  out  some  interesting  features.  In  particular,  the  present  results 
demonstrate  (for  the  clamped,  pressurized  spherical  cap  and  the  unrestrained 
spherical  cap  under  concentrated  load)  the  existence  of  a  number  of  bifurca¬ 
tion  points  in  the  axisymmetric  solution.  We  observe  that  although  two  dis¬ 
tinct  equilibrium  positions  are  possible  at  such  a  point,  the  distributions 
of  displacements  along  a  meridian  were  found  to  be  quite  dissimilar,  and  a 
change  from  one  path  to  the  other  would  involve  finite  (rather  than  infinitesi¬ 
mal)  changes  in  deflection.  It  is  interesting,  therefore,  to  note  that  the 
shell  under  concentrated  load  apparently  ignores  these  bifurcation  points 
since  the  load-deflection  curves  obtained  experimentally  by  Evan-Iwanowski 
et  al,12  are  in  very  good  agreement  with  the  branch  producing  the  first  maxi¬ 
mum  on  the  continuous  curve. 

In  the  case  of  uniform  pressure,  it  is  interesting  to  observe  that  if  one 
defines  the  critical  pressure  as  that  corresponding  to  the  first  maximum  on 
the  load-deflection  curve,  the  results  obtained  in  the  present  study  agree 
very  closely  with  the  four  most  widely  accepted  sets  of  results  for  the  axi- 

2-5 

symmetric  treatment  of  this  problem.  Furthermore,  our  results  for  the 
minimum  pressure  in  the  postbuckled  range  agree  quite  closely  with  those 
obtained  by  Thurston,5  who  was  able  to  obtain  results  in  this  zone  without 
establishing  the  continuity  of  the  global  solution.  One  clear  advantage  of 
having  the  complete  solution  is  the  following:  As  Budiansky2  points  out, 
referring  to  cases  where  the  maxima  were  not  located  with  certainty:  "Strictly 


11 


speaking,  it  must  be  conceded  that  the  upper  bounds  in  these  cases  have  not 
been  rigorously  established,  since  it  is  conceivable  that  failure  to  converge 
might  be  due  to  some  unknown  cause  other  than  the  nonexistence  of  an  adjacent 
equilibrium  position."  The  present  results  make  it  clear  that  the  results  of 
References  2  through  5  are  actually  the  first  maxima  on  the  load- deflect  ion 
curves.  They  also  reinforce  (though  any  further  reinforcement  is  hardly 
necessary  in  view  of  the  results  of  Reference  15)  Budiansky's  assertion  that 
one  must  Include  asymmetric  effects  in  the  study  of  the  buckling  of  spheri¬ 
cal  shells.  Ve  conjecture  here  that  the  axlsymmetric  global  solutions  mey 
prove  useful  in  the  more  difficult  problem  area  of  large  nonsymnetric  deforma¬ 
tions  by  providing  basic  states  about  which  to  expand,  particularly  in  the  post 
buckled  zone. 

It  is  considered  worth  mentioning  that  strain  and  potential  energies 
were  computed  for  the  various  deformation  states  and  the  so-called  energy 
criteria  of  buckling  (i.e.,  constant  volume  and  constant  pressure  buckling) 
were  applied.  The  results  again  agreed  quite  closely  with  those  reported  by 
Thurston,7  and  are  not  reported  here  in  the  interest  of  brevity.  However,  we 
observe,  strictly  speaking,  that  in  order  to  employ  such  an  energy  criterion 
one  should  have  the  complete  global  solution  available  in  order  that  the 
energy  levels  of  all  competing  equilibrium  states  be  compared. 

Finally,  we  remark  that  the  degree  of  compatibility  between  the  results 
of  the  complete  equations  and  the  small- finite  deflection  equations  is  some¬ 
what  surprising  since  the  latter  contain  only  the  simplest  type  of  nonlinearity 
whereas,  the  nonlinearity  in  the  complete  equations  is  more  complex.  The 
analytical  and  practical  advantages  offered  by  the  simpler  equations  are  suf¬ 
ficiently  large  -  for  example,  the  computer  program  for  the  simpler  set  takes 
only  one-third  as  long  to  run — that  exploitation  should  be  made  of  this  agree¬ 
ment  whenever  possible. 
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Figure  3c.  LOAD-DEFLECTION  CURVE  (UNRESTRAINED  EDOE:  K2 
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Figur*  4.  LOAD-DEFLECTION  CURVE  (CLAMPED  EDGE;  K2  •  144) 
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